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ABSTRACT 




The low energy spectra of protons and helium nuclei, recently measured 
by McDonald and Ludwig, and by Fan, Gloeckler and Simpson, are corrected 
herein for solar modulation and diffusive passage through interstellar matter. 
The resulting source spectrum is exhibited and compared with the spectrum 
of accelerated particles inferred from high energy data. This analysis yields a 
transmission function for the source environment which is similar in form to 
the velocity dependent solar modulation of low energy particles. 
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INTRODUCTION 

In the study of galactic cosmic radiation the evolution of the observed 
particles can be considered in four phases 1 : 

(a) injection from source regions 

(b) possible partial confinement in source regions 

(c) motion in the interstellar medium 

(d) solar modulation. 

Interstellar propagation has been studied extensively by many workers 2,3,4 
from a detailed determination of the charge spectrum of the cosmic radiation. 

In particular, it has been found that Li, Be, B (called L nuclei) occur in primary 
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cosmic radiation with very large abundance when compared to their universal 
abundance. As these nuclei are easily destroyed in nuclear interactions at 
stellar temperatures, it is reasonable to assume that they are probably not 
present in the source regions, but are produced in nuclear collisions of heavier 
nuclei with the protons of the interstellar medium. Using the abundance of 
these light nuclei relative to heavier nuclei, an estimate of the mean amount of 
matter traversed in the interstellar medium 5 has been made. For relativistic 
energies, a recent estimate 6 which takes into account effects due to decay of 
unstable isotopes formed in collisions and other details of nuclear collisions, 
gives a value of 2.5 gm/cm 2 as the mean amount of matter traversed by cosmic 
rays in the interstellar medium and is considered to be a reliable measure of 
this quantity. At low energies, however, the picture is not so clear and some 
experiments 8,9 seem to indicate an increase in the ratio of the abundance of L 
nuclei to that of heavier nuclei. Kaplon and Skadron 1 have interpreted the in- 
crease of this ratio as evidence that the amount of matter traversed by low 
energy nuclei is enhanced by a preferential confinement for low energy particles 
inside source regions that are surrounded by partially reflecting boundaries. 

The analysis presented herein is based primarily upon recent experiments, 
conducted on board the IMP-I satellite, which yielded the first clear measurement 
of the low energy spectrum of cosmic rays. The proton spectrum thus obtained 
by McDonald and Ludwig 10 and the helium nucleus spectrum obtained by Fan, 
Gloeckler and Simpson 11 are found to exhibit pertinent information about source 
transmission and interstellar propagation of low energy particles. Since these 
particles lie in the interval 20-100 MeV per nucleon, ionization loss is the 
dominant feature of their propagation in the interstellar medium, and their 
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energy is thereby greatly modified. This degradation of energy is a direct 
sensitive indication of the quantity of matter traversed by cosmic rays. Owing 
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to their low magnetic rigidity, these particles also serve as a probe of the 
magnetic structure that is characteristic of the source environment. 

In this paper we define the source spectrum in an operational way. We 
first take the experimentally observed energy spectrum and "demodulate" it 
for solar cycle effects. This "demodulated" spectrum which represents the 
cosmic ray energy spectrum outside the solar system, is traced back through 
the interstellar medium by taking into account the energy loss due to ionization. 
The path length dispersion arising from the diffusion of cosmic rays in the inter- 
stellar medium is considered in the analysis; this involves a specification of 
the distribution in path lengths about the mean value of 2.5 gm/cm 2 . Finally, 
by comparing the source spectrum obtained in this way with the spectrum of 
accelerated particles inferred from high energy data, we are able to reconstruct 
the transmission function for the physical source environment. 

THE LOW ENERGY PROTON AND HELIUM 
NUCLEUS SPECTRUM 

Figure 1 exhibits the proton and helium nucleus rigidity spectra obtained 
from recent balloon and IMP-I satellite measurements. The helium nucleus data 
have been normalized to the proton data using a scale factor of 7. Although the 
proton and helium nucleus rigidity spectra are comparable at high rigidities, it 
is clear that at rigidities below about 1500 MV the two spectra are definitely 
distinct. This shows that rigidity can not be used as a universal parameter for 
characterizing the spectrum of cosmic ray particles, at low energies. 
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Figure 2 shows the proton and helium nucleus spectra in a representation 
based upon the kinetic energy per nucleon. In this instance, the helium nucleus 
data have been scaled by 5.7 to achieve a normalization to the proton data. As 
indicated in Figure 2, a single curve describes all the pertinent proton and 
helium nucleus data points. This suggests that the kinetic energy per nucleon, 
or any other purely velocity dependent function, might be used to completely 
characterize the cosmic ray spectrum at low energies. The analytical expres- 
sion for the indicated curve is 

H = 10 8 E 1 - 5 (E + 500)- 4 (1) 

dE 

where dJ/dE has the units of particles/m 2 -sec-ster-MeV, and E is the kinetic 
energy per nucleon in MeV. This analytical representation of the observed 
spectrum is shown again as curve A in Figure 3. It was constructed to meet 
an asymptotic constraint imposed by the very high energy data 12 
(i.e. E > 10 10 eV); viz, 


« E- 2 - 5 . (2) 

dE 

The ratio of the proton flux to the helium nucleus flux at high energies is 
not as yet established with a precision comparable to that ascribed to this ratio 
for the low energies considered here. However, a comparison of the best 
estimates 13 * 14 of the cosmic ray flux at kinetic energies of 10 11 to 10 13 eV per 
nucleon with helium nucleon flux measurements 15 * 16 indicates that the ratio at 
very high energies is consistent with the value used in this paper for energies 
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in the interval 20 MeV per nucleon to 1 BeV per nucleon. In the region 2 to 15 
BeV per nucleon, integral intensity measurements give some evidence 17 that 
the proton to alpha ratio might be different at these intermediate energies, but 
this region is not crucial to the present analysis. 

For protons and helium nuclei in the low energy range the important ef- 
fects to be taken into account are ionization loss in the interstellar medium and 
solar modulation. At low energy these processes affect the protons and helium 
nuclei in an identical way as regards changes per nucleon. The reasons for this 
are that: (i) for equal Z 2 /A, which is the case for protons and helium nuclei, 
the ionization loss per nucleon is a universal function of velocity, and (ii) the 
low energy solar modulation is a universal function of velocity. 

McDonald and Ludwig 10 have deduced the galactic proton spectrum expected 
at Earth under the assumption of: 

(i) an energy spectrum at injection of the form 

/x(e) = 3" 77x 10" 4 ^8(1 + e ) 2 * 5 pro tons/m 3 -BeV, (3) 

where/x(e) is the density of cosmic rays of kinetic energy e in BeV and ft is the 
proton velocity in units of the velocity of light; 

(ii) 2.5 gm/cm 2 of hydrogen traversed in the interstellar medium; and 

(iii) a solar modulation given by 
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where (dj/ck) and (dj/de) 0 are the modulated and unmodulated spectra, respectively. 

The results of these calculations by McDonald and Ludwig are shown by the 
dotted curve of Figure 1. By comparing their proton data with this dotted curve, 
they conclude that either the source spectrum is much steeper than what has been 
assumed or the solar modulation is weaker at the low energies experimentally 
observed. Using the same solar modulation, the present analysis indicates that 
the source spectrum at low energies is strongly affected by the transmission 
characteristics of the source environment. 

SOLAR MODULATION 

Parker 18 has described the 11 -year solar modulation of cosmic rays in 
terms of a quasi- stationary solution to a Smoluchowski generalized diffusion 
equation 19 for the process of a charged particle diffusing among the magnetic 
kinks convected away from the sun by the plasma wind. For energies sufficiently 
low such that the Larmor radius (R) is small compared with the linear dimension 
(<l) of the magnetic kinks, the solution may be expressed as 

^in = eX P (- K //3)> < 5 ) 

where yu in is the cosmic ray density at the orbit of Earth, jj, out is the cosmic ray 
density outside the solar environment, and K is a solar modulation parameter 
that is proportional to the solar wind velocity and the thickness of the solar 
cavity, measured in units of kink intervals. 

For high enough energies such that R > Z the solution takes on the form 

= ($}■ 161 
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where Rg s Mc 2 /eB, M is the proton mass, B the magnetic field, and Z and A 
the particle charge and mass number respectively. 


Recent work on proton and helium nucleus spectra near solar minimum 
suggests that solar modulation effects are essentially velocity dependent. 20 The 
fact that protons and helium nuclei have the same low energy spectrum, in 
energy per nucleon, indicates that the dominant modulation is a function of 
velocity alone. For the present analysis, the change-over from the purely 
velocity dependent function (5) for low energies to the formula (6) for high 
energies is effected at a kinetic energy per nucleon (E„) given by 


E 0 = Me 2 



(7) 


Figure 3 shows the demodulated spectrum (curve B) thus obtained from the ob- 
served spectrum (curve A) for E Q = 1 BeV. Below approximately 500 MeV 
per nucleon, curve B exhibits the characteristic of a flat spectrum. We have 
found that curve B remains essentially the same for all (^Z/R 0 A) > 1. There- 
fore we may fit the data for helium nuclei (A = 2Z) as well as the data for 
protons (A = Z) by requiring that •£/ R 0 > 2. This may be achieved with reasonable 
values for £ and B (e.g. ~ 10 12 cm., B~ 10 -5 gauss). 


At intermediate energies of a few BeV the solar modulation could be a 
sensitive function of the detailed structure of the interplanetary magnetic field. 

In this paper, we rely primarily upon the purely velocity dependent solar modula- 
tion formula (5) for low energies and utilize the indicated hybrid formula (6) 
merely as a matching function that meets the asymptotic constraint of zero 
modulation at very high energies. 
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EXTRAPOLATION TO THE SOURCE 




By a detailed consideration of ionization losses, the solar demodulated 
spectrum is extrapolated to the source. Since these particles diffuse through 
the interstellar medium, the amount of matter traversed exhibits a dispersion 
about the mean value of 2.5 gm/cm 2 , and as these energies are low, the inferred 
source energy of each observed particle depends critically on the actual path 
traversed. An expression for the distribution of path lengths due to diffusion 
is derived in the appendix. The square root of the variance of the path lengths 
associated with this distribution is (2/5) 1/2 of the mean path length, and there- 
fore this complication to the extrapolation clearly cannot be ignored. 

Each observed particle is, in principle, mapped into a unit modulus energy 
distribution at the source according to the a priori distribution in path lengths 
arising from a diffusion process. These individual distributions at the source 
are weighted by the solar demodulation function, (5) and (6), and the resulting 
source spectrum (curve C) is shown in Figure 3. This source spectrum ex- 
hibits a peak at about 100 MeV per nucleon and falls rapidly at lower energies. 

DISCUSSION 

When compared with the observed spectrum, the source spectrum shows 
the systematic effects of energy loss by ionization, particularly at the low 
energy end. The smoothness at the low energy end is due to the dispersion in 
path lengths in the interstellar medium. If all particles were assumed to 
traverse exactly 2.5 gm/cm 2 , then the source spectrum thereby inferred would 
exhibit an abrupt cutoff at about 80 MeV per nucleon. 
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Under the assumption that high energy particles ( > several BeV) can 
freely propagate away from. the source accelerator and are not modulated 
drastically by galactic or solar magnetic fields, the energy spectrum of ac- 
celerated particles, at high energies, is the same as the asymptotic form (2) 
of the observed spectrum; 


£I=10 8 E- 2 * 5 . (8) 

dE 

This is shown as curve D in Figure 3. 

This spectrum (8) for accelerated particles probably does not remain 
valid down to extremely low energies, but we know that it is an adequate de- 
scription for several decades in the region of several BeV. As pointed out by 
Syrovatskii 21 this form (8) of the spectrum corresponds to a general situation 
of energy equipartition among cosmic rays, magnetic fields, and turbulent 
motion. If we assume that the same form of this spectrum is valid down to 
about 100 MeV, then the ratio of the source intensity (curve C) to the intensity 
given by curve D is a measure of the transmission efficency of the source en- 
vironment for the energy per nucleon under consideration. Figure 4 shows 
this transmission coefficient as a function of the reciprocal velocity. The fact 
that this coefficient exhibits a practically straight line in a semi-log plot with 
respect to 1//3 suggests that the accelerated source particles remain partially 
confined by a process of diffusion. The cosmic ray density "outside the source 1 
is thereby related to the density "inside the source" by a simple relation of 
the form 
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U(E)] out - k(E)l in exp (-K'//3) 


(9) 




where K' = 4.3, as calculated from the slope of the dashed straight line 
shown in the transmission ratio plot (Figure 4). 

Since the cosmic ray source spectrum appears to be a universal function 
of velocity rather than rigidity, we infer that the confinement to the source 
environment does not arise from trapping by ordered magnetic fields. Rather, 
a comparison of (9) with (5) indicates that the process involved here is one of 
diffusion, similar to that prevailing in the solar system. This suggests that 
the magnetic environment of the physical source of cosmic rays might be 
similar in form to the solar magnetic environment (e.g. magnetic turbulence 
convected outward by a plasma wind). 
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APPENDIX 


The dispersion in cosmic ray particle paths from the source to the solar 
system is here attributed to a diffusion mechanism. For a particle produced 
at the space -time origin of a homogeneous isotropic diffusive medium, the 
probability per unit time of first passage beyond a radius (a) is given by 14 

CD 

P(t) = 2 (^*) D exp ji(n + V)-n - Dtj , (Al) 

n=l V J 


where D is the diffusion constant characteristic of the medium. 


The mean time of first passage <t> is constructed as 


■/ 

''O 


<t> = tP(t)dt = a 2 /(6D). 


(A2) 


The second temporal moment <t 2 > is constructed as 


•'n 


<t 2 > = I t 2 P(t)dt = (7/180) (a 2 /D) 2 
'o 


(A3) 


The variance (<x t ) 2 of the first passage time is then given by 


(<r,) 2 . <(t -<t>) 2 > = <t 2 >-(<t» 2 = i- (a 2 /D) 2 


(A4) 
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By comparing the mean first passage time <t> , as given by (A2), with 


i 

the variance, as given by (A4), we note the relation 

<r t = (2/5) 1/2 <t> . { A5) 

If we assume that the quantity of matter traversed (S) by an observed 
particle is proportional to its total transit time, source to observer, then we 
can use the relation given by (A5) to infer that 

cr s s (<(S- <s>) 2 >) 1/2 = (2/5) 1/2 <s>- { A6) 

We consider a sample of observed cosmic ray particles which is char- 
acterized by a certain value for <( S {e.g. 2.5 gm/cm 2 ) and assign the cor- 
responding variance prescribed by diffusion theory. The well defined mean 
<( S y and square root of the variance a g (A6) here ascribed to (S) are sufficient 
to construct a normal Gaussian distribution function: 


P(S) 


( 277 )- 1/2 


exp 


(S - <s> ) 2 > 1 

2(a s ) 2 j 


(A7) 


Insertion of the explicit expression (A6) of the square root of the variance 
into (A7) yields 


P(S) 


( 5 / 77) 1 / 2 

2 <S> 


exp 



/ S - <s> 

V < s > 



(A8) 
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This distribution function (A8) is renormalized to account for the exclusion 
of negative values for (S), as follows: 


P' 



(A9) 


where 



Pds = 0.9429. 


This normalization procedure (A9) is based upon the approximation that 
is much greater than the shortest possible path length from the source to 
the observer. This assumption is valid in practice since the material diameter 
of the entire galaxy is an order of magnitude less than . 

The final expression of the renormalized distribution function is: 


P' (S) 


(5/n) 1/2 

— exp 

1.8858 <S> 



(A10) 


where 


I 


GO 


P' ds 


1 . 0000 . 
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FIGURE CAPTIONS 


Fig. 1 


Fig. 


Fig. 


Fig. 


Proton and He nucleus rigidity spectra. The He nucleus data have been 
multiplied by 7. The low rigidity data were obtained from the IMP-I 
satellite during the period November 1963— May 1964 (references 10, 
11). The other data were obtained from balloon measurements during 
June 1963 (references 22, 23). The dotted curve was calculated by 
McDonald and Ludwig (reference 10). 

The proton and He nucleus data (reference 10, 11, 22, 23) plotted in a 
kinetic energy per nucleon representation. The He nucleus data have 
been multiplied by 5.7. The solid line is equation (1). 

The energy per nucleon spectrum as observed (A), solar demodulated 
(B), corrected for ionization losses in interstellar matter (C), and 
extrapolated from higher energies, for comparison (D). 

Source transmission coefficient plotted as a function of 1 //3 (solid 
curve). The slope of the dashed straight line is -1.9. 
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